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D AUTO IONIZATION STATES OF THE HE AND H- 

A. K. Bhatia 
ABSTRACT 

Positions of the lowest 1,3 D e autoionization states of He and H" below the 
n = 2 level of the He + and H have been calculated va nationally, using Feshbach's 
Q-operator formalism. The trial wave function is of the Hylleraas-type with 
appropriate angular momentum factors. The widths and the shifts of the states 
have also been calculated. The shifts are found to be positive for all the states 
calculated here. The results with 112 terms for most states are lower than any 
previously calculated. The calculated lowest auto- ionization states ofthe He and 
H~ (relative to the ground states of He and H respectively) are 59.902 eV and 
10.1185 eV, in good agreement with the observed values of 59.9 eV and 10.13 ± 
0.015 eV. 


v 



D AUTIONIZATION STATES OF THE HE AND H 
I. Introduction 

A number of resonances 1 have been observed in He and H~ below the thre- 
shold of He + and H. Feshbach’s Q-operator formalism has been applied success- 
hilly to calculate 2, 3 the positions and the widths of the S and P autoionization 
states. The purpose of this paper is to extend these calculations to the D 
autoionization states in He and H" observed below the n = 2 threshold of the 
respective targets. These states lie in the continuum of electron scattering 
from single electron target system. They are not the stationary states of the 
two electron Hamiltonian and they autoionize by electron emission leaving behind 
bound states of the single electron system. The energy of the state can be 
written as 


E = S.Q+Aq (1) 

where £q is calculated variationally and Aq is the shift of £q due to the inter- 
action of the discrete state with the continuum. 

The most general D state wave function of even parity of two electrons 5 is 

$ (r lf r 2 ) = {(f ± f) [-J9f ( 0 , 6 0)+/ 3 cos J9f (0, fa 0 )] 

+ (f + ?) /3 sin & 12 J9^ - ( 6 , 0, 0)} (2) 

+ {(g ±g) [- cos 0 12 J9° + ( 6 , 0, 0 ) + /3 jQ ^ 0 ^]} 


1 



where the J9 are the rotational harmonics, depending on the symmetric Euler 
angles 6 ,4> , 4> 5 • These functions are eigenfunctions of exchange: indicating that 
they satisfy the following property 

8 12 ^ ±= ±( " (3) 

The trial wave function is of Hylleraas-type when the radial functions f = f(r 2 , 
r 2 , r ) and g = g (r lf r 2 , r 12 ) are given by 



It is implied in Eq. (2) that 


f = f C 2 . r,, r u ) (5a) 

and 

g = g ( r 2 - r i’ r i 2 > (5b) 

Using the properties (3) and (5) in Eq. (2), we see that the wave function is 
manifestly space-symmetric (upper sign) or space- antisymmetric (lower sign). 
The space symmetric and antisymmetric solutions correspond to singlet and 
triplet states respectively. 


2 



To best of our knowledge, this is the first calculation to use Hylleraas-type 
wave function for the D states. This calculation also shows the practical value 
of the symmetric Euler angle decomposition for higher angular momentum 
states. 

The eigenvectors in the two radial functions are independent. For simpli- 
fication in the calculation we have chosen 





The first term in Eq. (2) is formed by considering one electron in the s state 
and the other in the d state with total angular momentum equal to 2; the second 
term is formed by considering both the electrons in the p states with total 
angular momentum equal to 2. The other combination of the p and f states is 
included in the above wavefunction. The upper sign corresponds to the singlet 
states and the lower sign corresponds to the triplet states. 

The projection operator is given by 

Q = 1 - Pj - P 2 + Pj P 2 (6a) 

= 1 - P (6b) 

P i = 1 K ( 7 i> > <( t>o ( 7 i)l (7) 


3 



where <fi 0 (r . ) is the ground state of the ith particle in the nuclear field of 
charge z: 


<f!) (8) 

R„ (r) = 2z«e-" (9) 


The Q-operator given here is restricted to the autoionization states below the 
n = 2 threshold. The projection equality, Q 2 = Q is valid here. 

II. Calculations and Results 

The expectation value of the energy in the restricted Q<I> space is given by 
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Q- 


(OQHQ$) 
<<D Q <fr> 


( 10 ) 


The decomposition of Eq. (10) is discussed in the Ref. 2. In Table I we present 
results for the singlet D states of the He and H as a function of the Perkeris 
numbers for D states i.e., the number of terms 


CO 

N (co) = 2 y | n (co. ; 


( 11 ) 


co . = 0 


where n (ox) contains all terms r^ r^r^ 2 such that^ + m + n = 0, 1, 2, . . . 5. 
The factor 2 in Eq. (11) arises due to two types of terms (4a) and (4b) in 
Eq. (2). The nonlinear parameter y is kept fixed at 0.70 and 0.50 for He 


4 



and H states respectively. The other non-linear parameter S is varied to get 
the optimum Sq. Only one singlet state of H - has been calculated. The higher 
singlet and triplet states do not descend below the n = 2 threshold in this calcu- 
lation due to their enormous size. A general configuration interaction wave- 
function instead of Eq. (2) may be better for these states. In Table II we 
present results for the triplet D states. Some of the results for N = 8 and 20 for 
higher states in Tables I and II lie above threshold of n = 2, because they are not 
optimized for these terms. 

An exact formula for the shift A q of an isolated resonance can be written 
as follows 6 : 


a Q 


J_ f, f r(E') dE- 
2 77 I e - E' 


( 12 ) 


where <p indicates the principal value. The width of the resonance is given by 


r (E) = 2 k I <P T (E) I H | Q 3> > | 2 ( 13 ) 

In these formulas rydberg units are used throughout. Q3>is the exact eigenfunction 
of the projected problem and E is the total energy of the resonance. If the Hamil- 
tonian is written in the form 


H = Hq + 2/r 12 


(14) 


5 



with the observation that P and Q commute with H 0 and that PQ = 0,F (E) reduces 
to the form 


r (E) = 2 k 


pT (E)| 



(15) 


PT(E) is the solution of the optical potential problem less the resonant term 4 : 


(H' - E) P T i 0 (16) 

H' = H pp + U' (17) 


H pp is the exchange approximation Hamiltonian andl)' is the optical potentials 
less the resonant term. Various approximations 3 to the potential are possible 
Since the most important correction to the exchange- approximation phase shifts 
comes from the polarization 7 of the target, the continuum wave function T has 
been calculated by the method of polarized orbitals 7 : 


(poi) u (r ) (poi) 

T (E) (fij) [<£ 0 (r 2 ) + 4> (r v r 2 )) ■ 

r i 


(18) 


±( 1 - 2 ) 


The explicit inclusion of the z-factor for positive ion targets in the method of 
polarized orbitals is given in Ref. 8, where 


6 



e ( r r r 2> 


— zr 


(19) 


(p°£) 

^( r ,. r 2 ) 


2 (l 2 V° S ^12 

\2 z r2 + 7 sn 


and e (r y , r 2 ) is a step function. 

The function u (r) /r is normalized as a plane wave or its Coulomb counter- 
part: 


1 im u (r) = _ sin 
r -> 00 k 


r + cr 


4 




(20) 


where is the Coulomb phase factor: 




irg r ^ + 1 - 1 (2 k l ^j + 1—1 In (2 k r) z>l 


a l = ] 


( 21 ) 


z = 1 


The quantity 77^ is the residual phase shift. The energy of the scattered particle 
is k 2 and is related to the total energy E by 

E = E t + k 2 ( 22 ) 

where E t is the ground state energy of the target system (He + or H). 

In Tables III and IV we present our results for the positions, widths and 
shifts for singlet and triplet states respectively, and compare them with other 
calculations and experimental results. The positions of the resonance are given 
relative to the ground state of the target and the results include the reduced 


7 



mass correction,, The reduced rydberg is 13.60350 eV and 13.59794 10 eV for 
He and H - respectively. The comparison indicates that the present calculation 
gives positions of states lower than other calculations except for the 7th state in 
the singlet and triplet series for He. The shifts for all the states are found to be 
always positive. The shifts are always small and comparable to the widths. But 
in few triplet states of He, the shifts are much larger than widths, still 
small. The calculated lowest singlet autoionization state of He is 59.902 eV in 
agreement with the observed values of 59.95 11 eV and 59. 9 12 eV. The calculated 
lowest autoionization state of H is 10.1185 eV, in good agreement with the ex- 
perimental value 10.13 ± 0.015 eV of McGowan 13 et al. and is within the ex- 
perimental error. There are no experimental results available for higher singlet 
and triplet states. The widths of lowest singlet states of Ke and H~ agree with 
the results of Cooper 14 et al. and Burke and Taylor 1 respectively. 

In conclusion we will urge measurements of positions and widths of these 
states for comparison purposes. 
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Table I 


Auto- Ionization Energies (Ry) as a Function of N for 1 D States of He and H' 


N ( o >) 

He (y = 0 . 7 ) 

H " ( 7 = 0 . 5 ) 

-e, 

4 = 1.0 

-e* 

6 = 0.5 

- ? 3 

» = 0.4 

-8 

8= 6.4 

8 = 0.3 

" e . 

8 = 0.2 

-e, 

8 = 0.2 

- g 8 
8 = 0.2 

- e . 

8 = 0.2 

8 ( 1 ) 

1.403467 

1.093372 

1.068919 

0.943733 

0.910025 

0.652551 

0.502781 

0.002152 

0.255038 

20 ( 2 ) 

1.405219 

1.131226 

1.106968 

1.033154 

1.028373 

0.998156 

0.967890 

0.917407 

0.255848 

40 ( 3 ) 

1.405558 

1.137801 

1.112074 

1.062206 

1.057408 

1.053822 

1.024525 

1.009070 

0.256129 

70 ( 4 ) 

1.405622 

1.138604 

1.112777 

1.071807 

1.062111 

1.058384 

1.038589 

1.035262 

0.256160 

112 ( 5 ) 

1.405634 

1.138752 

1.112855 

1.073241 

1.062920 

1.058566 

1.043634 

1.039609 

0.256174 


Table II 


Auto- Ionization Energies (Ry) as a Function of N for 3 D States of He 



He (7= 0.7) 

N(<w) 

-£i 

8 = 0.4 

- g 2 

8 = 0.4 

- £ 3 

b — 0*3 

8= 0.3 

- £ s 

8= 0.2 

- e 6 

8 = 0.3 

_8 7 
8= 0.2 

- g 8 

8= 0.2 

8 (1) 

1.155374 

1.094690 

1.016989 

0.984820 

0.952665 

0.570581 

0.494916 

0.017835 

20 (2) 

1.165450 

1.117471 

1.069293 

1.043720 

1.033629 

0.966562 

0.968083 

0.917877 

40 (3) 

1.167327 

1.120846 

1.081119 

1.062914 

1.056817 

1.026804 

1.023591 

1.009192 

70 (4) 

1.167581 

1.121296 

1.083075 

1.066306 

1.058428 

1.046344 

1.038565 

1.034619 

112 (5) 

1.167611 

1.121361 

1.083336 

1.066831 

1.058605 

1.049449 

1.041638 

1.037735 













Table in 


Compilation of Theoretical and Experimental Results. Units: eV 


System 
and State 

Resonance 

Parameters 

This 

Calculation 

Cooper 
et al.° 

Altick 

and 

Moo re b 

Burke 
et al. c 

Burke 

and 

Taylor 3 

Ormonde 
et al." 

Expt. 

He *D(1) 

e 

59.8801* 


60.115 






r 

0.0729 

0.0732 

0.0748 


0.0662 




A 

0.0220 






(60.0 f ’ e 


E 

59.902 

60.025 



59.911 


ij 59.95 h 

* D (2) 

e 

63.5106 


63.601 




| S 9.9 1 


r 

0.0187 

0.0165 

0.0179 






A 

0.0044 








E 

63.515 

63.575 






*D(3) 

e 

63.8629 


63.904 






r 

5.807 x 10' 4 

2.83 X 10- 4 

4.04 X 10' 4 






A 

1.815 X 10’ 4 








E 

63.863 

63.897 






‘d<4) 

e 

64.4018 


64.480 






r 

7.124 X 10‘ 3 

7.1 x 10- 3 

11.6 x 10‘ 3 






A 

1.667 X 10' 3 








E 

64.403 

64.429 






*D(5) 

e 

64.5422 


64.638 






r 

4.045 X 10‘ 4 

1.68 x 10~ 4 

3.11 X 10- 4 






A 

1.045 X 10- 4 








E 

64.542 

64.557 


64.580 




*D(6) 

E 

64.6014 


64.682 






r 

5.104 X 10~ 7 








A 

1.237 X 10' 7 








E 

64.601 

64.611 


64.539 




D(7) 

e 

64.8045 








r 

1.777 X 10- 3 

3.68 X 10‘ 3 







A 

0.691 X 10' 3 








E 

64.805 

64.797 


64.796 




'D(8) 

e 

.64.8593 








r 

3.891 X 10" 4 

0.973 X 10- 4 







A 

1.090 x 10- 4 








E 

64.859 

64.861 


64.880 




H _ *D(1) 

e 

10.1145 








r 

0.0100 




0.0088 

0.0088 



A 

0.0040 








E 

10.1185 




10.125 

10.126 

10.130 ±.015' 


b 

c 

d 


Cooper, Ormonde, Humphrey and Burke, Ref. 14. 

P. L. A! tick and E. N. Moore, Proc. Phys. Soc., 92, 853 (1967). 

P. G. Burke, D. D. MeVicar and K. Smith, Phys. Letters 12, 215 (1 964) - 
P. G. Burke, Ref. 1. 


®$. Ormonde, J. McEwen and J. W. McGowan, Phys. Rev. Letters, 22, 1165 (1969). 

*M. Rudd, Phys. Rev. Letters J3, 503 (1964); J5, 580, (1965); J. A. Simpson, S. R. Mielciorek and J. W. Cooper, J. Optical Soc. 
Amer. 54, 269 (1964). 

®J. J. Quern dne'r, C, Paquet and P. Marmot, Phys. Rev. A4, 496 (1971). 

^ Burrow, Ref. 11. 

'Oda, Nishimura and Tahira, Ref. 12. 

•McGowan, Clarke and Curley, Ref. 13. 





























Table IV 


Compilation of Theoretical Results. Units: eV 


System 

and 

State 

Resonance 

Parameters 

This 

Calculation 

Cooper 
et ad. 8 

Altick 

and 

Moore 

Burke 
et al. c 

He 3 D(1) 

e 

63.1180 


63.157 



r 

2.715 x 10- 6 

1.44 x 10' 6 

0.12 x 10 -6 



A 

0.260 x 10- 3 





E 

63.118 

63.141 



3 D(2) 

6 

63.7472 


63.797 



r 

1.919 x 10- 4 

2.48 x 10" 4 

2.53 x 10- 4 



A 

0.917 x 10" 4 





E 

63.747 

63.796 



3D(3) 

e 

64.2645 


64.284 



r 

3.305 x 10- 6 

0.468 x 10- 6 

5.3 x 10' 6 



A 

1.203 x 10 " 4 





E 

64.265 

64.273 



3 D(4) 

e 

64.4890 


64.559 



r 

1.362 x 10' 4 

0.128 x 10' 4 

2.05 x 10- 4 



A 

0.502 x 10~ 4 





E 

64.489 

64.509 


64.486 

3 D(5) 

e 

64.6009 





r 

1.192 x 10- 7 

0.79 x 10- 7 




A 

5.101 X 10~ 7 





E 

64.601 

64.610 

64.665 

64.562 

*0(6) 

e 

64.7255 





r 

8.899 x 10- 8 





A 

6.656 x 10~ 5 





E 

64.726 



64.710 

*0(7) 

e 

64.8317 





n 

l 

5.967 x 10" 5 

7.32 x 10" 5 




A 

1.173 x 10- 5 





E 

64.832 

64.836 


64.823 

3 D(8) 

e 

64.8852 





r 

2.596 x 10- 7 





A 

1.43 x 10 - 7 





E 

64.885 





“Cooper, Ormonde, Humphrey and Burke, Ref. 14. 

b P. L. Altick and E. N. Moore, Proc. Phys. Soc.,^2, 853 (1967). 

C P. G. Burke, D. C. McVicar and K. Smith, Phys. Letters J2, 215 (1964). 
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